 **********************************
 2D FEM v54.3 of 22/08/18
 **********************************
 Ne =           20
 Nn =           42
 Nb =           42
 Nm =            1
 Np =           20
 Nd =            4
 Ng =            2
 Ns =            3
 ldw =          336
 Sx, Kappa, tc, por =    1.00000000000000       0.500000000000000     
   1.00000000000000        1.00000000000000     
 tmax, dt, dto, Nt
   1500.00000000000       0.100000000000000        300.000000000000     
       15000
 theta=    1.00000000000000     
 ipar =            0           0           0         336           0          -1
           0           0           0           0           0           0
           0           0           0           0
 fpar =   1.000000000000000E-009  1.000000000000000E-012  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000
 maxit, Ctol =        1000  1.000000000000000E-003
 i, xg(i), wg(i)
           1 -0.577350269189626        1.00000000000000     
           2  0.577350269189626        1.00000000000000     
 maximum advection   CFL (Vabs * dt / h)       =   4.753328917976242E-004
 maximum dispersion  CFL (Dabs * dt / h**2)    =   3.819249589742163E-005
 maximum Peclet number (0.5 * Vabs * h / Dabs) =    6.22285714285715     
 SUPG parameter, aopt =   0.839309981158420     
 Initialization time=  1.5600100E-02  seconds
 Total time=   1.014007      seconds
